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A promising technique for measuring single electron spins is magnetic resonance force microscopy 
(MRFM) , in which a microcantilever with a permanent magnetic tip is resonantly driven by a single 
oscillating spin. If the quality factor of the cantilever is high enough, this signal will be amplified over 
time to the point that it can be detected by optical or other techniques. An important requirement, 
however, is that this measurement process occur on a time scale short compared to any noise 
which disturbs the orientation of the measured spin. We describe a model of spin noise for the 
MRFM system, and show how this noise is transformed to become time-dependent in going to the 
usual rotating frame. We simplify the description of the cantilever-spin system by approximating 
the cantilever wavefunction as a Gaussian wavepacket, and show that the resulting approximation 
closely matches the full quantum behavior. We then examine the problem of detecting the signal 
for a cantilever with thermal noise and spin with spin noise, deriving a condition for this to be a 
useful measurement. 

PACS numbers: 68.37.Rt, 33.35.+r, 07.79.Pk 



I. INTRODUCTION 

Sidles, in 1991, first proposed the use of magnetic reso- 
nance and mechanical oscillators to sense a weak force [l| . 
Ever since this proposal, there has been considerable 
progress in the experimental implementation of such a 
technique culminating in single spin detection in 

2004 Q. That the force detected is of the order of at- 
tonewtons, highlights the usefulness of magnetic reso- 
nance force microscopy (MRFM) based techniques. 

MRFM based techniques have a wide range of applica- 
tions. One such application is imaging; in fact, imaging 
was the original motivation for proposing MRFM; imag- 
ing at the nanoscale like that of biological molecules such 
as proteins, viruses etc., is of immense value to the soci- 
ety at large [E^. Another application, and the one that 
is of interest to us, is the use of MRFM for single spin 
measurements. Numerous models of quantum computing 
involving spins have been proposed (6|, and these mod- 
els require the ability to do single spin measurements, in 
addition to purely scientific value of such a measurement 
capability. 
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In 2001, Stipe et al. [s^ first implemented the OScil- 
lating Cantilever-driven Adiabatic Reversals (OSCAR) 
protocol. The frequency shift that takes place in this 
protocol, can be measured experimentally with high pre- 
cision. The idea of this spin manipulation protocol is to 
transform the cantilever-spin interaction force into a shift 
in the resonant frequency of the oscillating cantilever, by 
using a gain-controlled feedback mechanism; the inter- 
action force between the cantilever and the spin, which 
is either attractive or repulsive depending on the orien- 
tation of the spin, gets transformed to a positive or a 
negative shift in frequency; by measuring this shift one 
can determine the orientation of the spin. 

Theoretical analysis quickly followed with a series of 
work by Berman et al. [7-10] and Brun and Goan [lL[l2j. 
Berman et al. use a semi-classical approach in their 

analysis, wherein they treat cantilever as a classical Har- 
monic oscillator while treating the spin quantum mechan- 
ically. Brun and Goan [III, [HI provide a fully quantum 
description of the system, in that they treat both can- 
tilever and spin quantum mechanically; they also con- 
sider decoherence effects acting on the cantilever in their 
analysis, including continuous measurement of the can- 
tilever due to optical inter ferometric techniques. Though 
modeling and analysis of the system using a fully quan- 
tum description is an important theoretical development, 
Brun and Goan [III, [Hf describe the system as a pure 
quantum state, instead of a (more general) density ma- 
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trix, and use quantum (pure) state diffusion techniques 
to describe their evolution, instead of a stochastic mas- 
ter equation evolution [13|; also, they do not consider 
decoherence effects acting on the spin. 

In this work, we analyze the MRFM system in gener- 
ality within a Markovian framework. We describe the 
system using a density matrix, and evolve it using a 
stochastic master equation [13]. We consider decoher- 
ence effects acting on both the cantilever and the spin. 
The stochastic master equation, though it gives a fully 
quantum-mechanical description of the system, is numer- 
ically expensive. To ameliorate this, we show that the 
state of the cantilever can be approximated as a Gaus- 
sian wave packet, leading to a description of the MRFM 
system by a closed set of 11 coupled stochastic differential 
equations. We show that the Gaussian approximation is 
valid in the parameter regime of interest, by numerically 
comparing the evolution of the Gaussian equations with 
the fully quantum stochastic master equation. 

Further, we use the Gaussian equations to analyze the 
OSCAR protocol as a tool for single-spin measurement. 
We consider the constraints set by the spin noise on a 
single-spin measurement. For the parameter values cho- 
sen, we calculate a bound on the spin noise time scale. 
We choose a rate for spin noise that satisfies this time 
scale, and numerically show that the OSCAR protocol 
can indeed be used to do single-spin measurement. 



A. Overview of the paper 

Sec. |II] gives an overview of the MRFM system and 
the OSCAR protocol. In Sec. Hill we use the adiabatic 
approximation and derive an effective Hamiltonian for 
the OSCAR MRFM system. We consider the decoher- 
ence effects acting on both the cantilever and the spin, 
in Sec. |IVl the sources of decoherence include the ther- 
mal bath, continuous measurement of the cantilever us- 
ing optical interferometry, and spin noise due to mag- 
netic sources. We then discuss the feedback mechanism 
used to implement the OSCAR protocol (Sec. |V]). We 
derive equations for the moments in Sec. IVH we consider 
the unitary case first (Sec. IVI Ap . before proceeding to 
the more general open system evolution (Sec. IVIBp . Fi- 
nally, we present the numerical results in Sec. IVIII We 
first discuss the parameter values that are used in our 
simulations; we consider the effectiveness of the Gaus- 
sian approximation by comparing it to the fully quan- 
tum stochastic master equation evolution (Sec. I VII A^ : 
and lastly, we ask what it takes for an OSCAR MRFM 
system to be an useful single-spin measurement device. 



II. OSCAR MRFM - MODEL AND OPERATION 

The basic model of MRFM involves a cantilever oscil- 
lating in close proximity to a spin that is contained in 
a substrate (Fig. [1]). A ferromagnet that is attached at 



the tip of the cantilever interacts with the spin. This 
interaction changes the amplitude and frequency of the 
oscillator, which is measured to determine the orientation 
of the spin. 

There are three different sources of magnetic field 
present in the system: (i) a constant field, Bq^ in the 
z direction, (n) an oscillating ac microwave (RF) field, 
Bac^ that is applied in the x — y plane, and, {Hi) a dipole 
field due to the ferromagnetic tip, B^. 

The constant field Bq polarizes the spin in the z direc- 
tion. The rotating microwave (RF) field, together with 
the dipole field from the ferromagnet, causes the spin 
precession axis to periodically reverse with the same fre- 
quency as the cantilever. When in resonance, the spin os- 
cillates with the same frequency as that of the cantilever. 
Under the resonant condition, the force due to interac- 
tion of the ferromagnet at the tip of the cantilever and 
that with the magnetic moment of the spin gets ampli- 
fied: this amplification leads to a change in the resonant 
amplitude and frequency of the cantilever. This change 
can be observed by monitoring the cantilever, usually 
through optical interference. 

In the OSCAR protocol, the oscillation of the can- 
tilever (ferromagnetic tip) causes adiabatic reversals of 
the spin, which in turn interacts with the cantilever to 
change the resonant amplitude and frequency of the can- 
tilever Q. In this protocol, feedback is used to shift the 
resonant frequency of the cantilever; the change in res- 
onant frequency — positive or negative — of the cantilever 
is used to determine the orientation of the spin. The 
feedback mechanism used is positive gain controlled, and 
is used to maintain the cantilever amplitude at a pre- 
determined fixed value. The cantilever thus behaves as a 
frequency determining element p^ . 



III. EFFECTIVE HAMILTONIAN USING THE 
ADIABATIC APPROXIMATION 

The cantilever and spin form the two parts of the sys- 
tem, and are described by their respective Hamiltonians: 

Hc{t) = film + mcj2zV2 - f(t)Z 

Hs{t) = eS, - r]ZS, ^ ^ 

where He if) and Hs{t) are the Hamiltonians of the can- 
tilever and spin, respectively [9]; f{t) is the positive gain 
controlled feedback mechanism that implements the OS- 
CAR protocol (refer Sec.[V]); and, e and r] represent the 
strength of the microwave (RF) field and spin cantilever 
interaction, respectively. Note that we have included the 
interaction term rj Z Sz as part of the spin Hamiltonian 
Hs{t). 

In the parameter regime of the OSCAR protocol, the 
time scales of the spin and the cantilever are well sepa- 
rated. The frequency of spin precession is much greater 
than the oscillation frequency of the cantilever. On the 
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FIG. 1. (Color online.) MRFM system: model and operation. A cantilever with a ferromagnetic tip oscillates in close proximity 
to a substrate with a free spin. A uniform magnetic field Bq is applied in the z direction and a microwave (RF) field is applied in 
the x-y plane. The position of the cantilever is monitored continuously by shining a laser on the cantilever tip and using optical 
interference to make measurements. A feedback mechanism (not shown in the figure) maintains the amplitude of the cantilever 
at a pre-determined fixed value, making the cantilever a frequency-determining element. The direction of the frequency shift 
reveals the orientation of the spin. 



time scale of the spin, the cantilever part in the interac- 
tion term in the Hamiltonian, 

Hi = -vZS, 

can be treated as a constant; thus, the quantum can- 
tilever position term, Z, on the time scale of the spin can 
instantaneously, in the adiabatic limit, be treated as a 
classical function Z. 

The time dependent spin Hamiltonian in the adiabatic 
limit is given by 



Hs{t) = eS, - TjZS,, 
whose instantaneous eigenstates are 



(2) 



\v±{t)) 



+ (±A + r?Z)2 

±X + riZ 



I i>, (3) 



+ (±A + 7?Z) 
with eigenvalues 

X{t) = x/e" + r?2z2. (4) 

If the initial state of the spin is 

IVs(O)) = a(0)|t;+(0)) + /?(0) |t;_(0)), (5) 

then the state of the spin at later times can be written 
as 



\Mt)) = a{t)\v+{t)) + m\v-{t)), 



(6) 



where \v±{t)) is the instantaneous eigenstate of the time- 
dependent spin Hamiltonian at time t. The time evolu- 
tion of the state is given by the Schrodinger equation 



d_ 

di 



Hsit)\Mt))- 



(7) 



Using Eq. (|6|) in Eq. ([7]), and solving for a{t) and /3{t) 
using the adiabatic approximation, we get 

a{t) = a(0)exp|--^^ X{s) ds^ , (8) 

m = /3(0)exp|^£ A(8)d5|. (9) 

The probabilities of the spin being up or down in the 
reference frame of the effective magnetic field are con- 
stants of motion, and are given by |q^(0)P and |/3(0)p, 
respectively. In the solution above, we have used the fact 
that the time scale of cantilever motion is much slower 
than the time scale of spin precession, i.e., the rate of 
change of the cantilever \dZ/dt\ <C (e^/r^), and conse- 
quently, we have ignored terms proportional to dZ/dt in 
deriving Eqs. (j8j) and ([9]). 

Assuming the cantilever and spin are initially decou- 
pled, the initial state of the cantilever-spin system is 



(10) 



where \ipc{^)) is the initial state of the cantilever, and 
ItpsiO)) is the initial state of the spin as defined in Eq. (j5j). 
Further, we assume that the cantilever is localized, and 
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consequently the motion of the cantilever can be de- 
scribed by two wave packets — one corresponding to each 
spin orientation. 

As we know, \v±{t)) are the eigenstates of the spin 
Hamiltonian, and satisfy 

Hs{t)\v±{t)) = ±X{t)\v±{t)), 

where X{t) is its eigenvalue [Eq. (j4])]. We assume that 
the microwave (RF) field, e, is much stronger than the 
interaction strength between the cantilever and the spin 
rjZ. Expanding A(t) as a binomial series, we get 

X{t) = 6 [1 + {vZ/ef/2 + 0{{7iZ/ef)] . (11) 

As the microwave (RF) field dominates the interac- 
tion strength, we can ignore third — and higher order — 
contributions in this expansion. Since the first term in 
the expansion is a constant, it amounts to a phase factor, 
and can be ignored as well. 

Let us now consider the action of the total Hamiltonian 
on the joint cantilever and spin states: 

{Hc + Hs{t)){\z)(E)\v±{t))) = {Hc±m){\z)®\v±{t))), 

(12) 

where 1^) is the position of the cantilever at time t. We 
now quantize the position function to obtain the effective 
Hamiltonian acting on the cantilever: 

^' _ { He + FZ^, if spin-up, i.e. 
^ ~ \ Hc - FZ^, if spin-down, i.e. \v-{t)) 

= He ^ 2T a'z (13) 

where we have set 

F = r/V2e, (14) 

and (j^ is the rotating operator that represents the spin 
with respect to a reference frame that moves with the 
effective magnetic field. The total effective Hamiltonian 
acting on the cantilever + spin system is, thus, 

H = He- (15) 

IV. MODELING DECOHERENCE IN OSCAR 
MRFM SYSTEM 

In reality, quantum systems are never completely iso- 
lated. The interaction of the system with the rest of the 
universe is seen as decoherence on the system of interest. 
There are various decoherence processes that act on both 
the cantilever and the spin part of the system. 

A. Decoherence in the cantilever 

Two primary sources of decoherence act on the can- 
tilever: (a) the thermal environment, and, (6) the con- 
tinuous measurement of the external read out. Brun and 



Goan [ll| have shown that these processes can be mod- 
eled as Lindblad operators of the form, 

Lk = AkZ ^ iBkP, (16) 

where Li captures the effect of thermal environment and 
L2 is due to continuous monitoring of the cantilever po- 
sition. The coefficients are 

Ai = ^A -jrnmkBTj h? , Bi = ^-irn/^rnkBT 
A2 = ^^K^Whl, B2 = 0. 

(17) 

Modeling the thermal environment as a Lindblad opera- 
tor Li requires an additional damping term to be added 
to the Hamiltonian (refer [ll| Sec. IV), and the new ef- 
fective Hamiltonian of the cantilever becomes 

Heff = H ^ ImR. (18) 

where H is the effective Hamiltonian from Eq. (p!5]) . 7m 
is the rate of damping due to the thermal environment, 
and 

R = ^{Zp + pZ). (19) 

B. Modeling spin noise 

The source of noise acting on the spin could either be 
dipole-dipole interactions of the spin with other spins in 
the lattice, or noise in the magnetic field. The dipole- 
dipole interaction acts in the laboratory frame while the 
noise due to magnetic sources can be treated in the ref- 
erence frame of the effective magnetic field. 

Thermal noise is one of the dominant decoherence pro- 
cess in the MRFM system. This noise affects the can- 
tilever motion, which in turn affects the interaction be- 
tween the cantilever and the spin. Thus, noise in the 
cantilever motion means that there is noise in the dipole 
magnetic field due to cantilever-spin interaction. Also, 
the external microwave (RF) field is not free of noise. 
Stipe et al. [3| have shown that the spin-lattice relax- 
ation time is quite long, on the order of seconds. This 
means that spin decoherence is dominated by magnetic 
noise. As a result, we consider spin relaxation due to 
magnetic noise and ignore spin-lattice relaxation. 

In the reference frame of the effective magnetic field, 
spin relaxation due to magnetic noise corresponds to sim- 
ple spin-fiip noise. The Lindblad operator corresponding 
to spin noise can thus be modeled as, 

L3 ^/K^ax' (20) 

Kg is the rate of spin noise due to magnetic sources, and 

4 = \v+{t)){v.{t)\ + \v.{t)){v+{t)l (21) 

where \v±{t)) are the instantaneous eigenstates of the 
spin Hamiltonian at time t [Eq. (|3|)]; note that spin- up 
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or spin-down here corresponds to spin in the direction of, 
or opposite to, the effective magnetic field. 

For MRFM to be useful as a single spin measurement, 
it is necessary that the spin decoherence rate tZg be small 
compared to other parameters in the cantilever-spin sys- 
tem. As we measure the cantilever, we learn something 
about the state of the spin, and eventually, the spin re- 
laxes to one of its eigenstates, either in the direction of 
the effective field or in the opposite direction. There are 
two different time scales at work here: (i) the spin lo- 
calization time, and, (ii) our observation time; the spin 
localization time scale is the time it takes, on average, for 
the spin to collapse to one of its eigenstates due to con- 
tinuous monitoring of the cantilever; on the other hand, 
the observation time scale is the time it takes for us to 
'know' the state of the spin. Typically, the observation 
time scale is greater than the spin localization time scale 
as it would take a little longer after the spin is completely 
relaxed to 'know' the state of the spin; however, one can 
also conclude the state of the spin based on the trend 
(depending on SNR) in the frequency shift. In the latter 
case, we do not wait for the spin to localize completely 
to guess the state of the spin. Irrespective of the method 
used, the time scale of spin noise has to be longer 

than both spin localization and observation time scales 
for MRFM to be useful in single-spin measurement. 



the cantilever amplitude at a pre-determined fixed value. 
The feedback has the form 

f{t) = g X {AMP - Amp{t)) x I,{t - A), (23) 

where g is the feedback gain, AMP is the pre-determined 
set point amplitude, Amp{t) is the cantilever amplitude 
at time t, Ic{t — A) is the delayed output photocurrent, 
and A corresponds to a delay of 7r/2 radians, or equiv- 
alently, a fourth of the cantilver oscillation time period. 
The current amplitude of the cantilever, Amp{t)^ is de- 
rived from the measured output photocurrent Ic{t)^ using 
simple signal processing techniques. 



VI. MOMENT EQUATIONS FOR OSCAR 
MRFM WITH GAUSSIAN APPROXIMATION 

In this section, we derive moment equations for the 
cantilever and spin. We consider the unitary case first 
and then move to open system evolution. We use the 
Schrodinger equation to describe the unitary evolution, 
and a full quantum stochastic master equation for the 
open system evolution, with continuous monitoring of the 
cantilever motion. Finally, we make Gaussian approxi- 
mation for the cantilever degree of freedom and show that 
the moments then obey a closed set of coupled equations. 



V. CONTINUOUS MEASUREMENT AND 
FEEDBACK IN OSCAR PROTOCOL 



A. Unitary evolution 



In the OSCAR protocol, the spin orientation is mea- 
sured by measuring the frequency shift of the cantilever. 
A positive gain-controlled feedback mechanism maintains 
the amplitude of the cantilever at a pre-determined con- 
stant, leading to a change in cantilever frequency. 

The continuous monitoring of the cantilever motion 
is done by optical interferometry. As shown in Fig. (TJ 
a laser is placed close to the tip of the cantilever; an 
optical microcavity is formed with the cantilever on one 
side and the cleaved end of the fiber (laser) on the other. 
Since the motion of the cantiever is slow compared to 
the optical frequency, this system can be analyzed in the 
adiabatic limit. A homodyne measurement is carried out 
on the light that escapes this cavity. The output of the 
homodyne measurement corresponds to the position of 
the cantilever (Z), and is given by 



Icit) = /? 



SedKE^ dWt 

{Z} + y/jced —7- 

7c dt 



(22) 



Here, Ic{t) is the output photocurrent of the homodyne 
measurement, n is the coupling between the cantilever 
and the cavity, 7c is the cavity loss rate and ed is the 
detector efficiency. Brun and Goan [H] (Sec. V) have 
analyzed this system in detail and we refer to that paper 
for more information. 

The feedback mechanism in the OSCAR protocol is 
positive gain-controlled, and its objective is to maintain 



The evolution of the cantilever-spin system is given by 
a Schrodinger (or a von Neumann) equation: 



dp 
~dt 



(24) 



H is the effective Hamiltonian of the cantilever-spin sys- 
tem [Eq. (p^ ]. 

To derive the moment equation for an arbitrary op- 
erator O, we apply the Schrodinger equation to obtain 
equation for (O): 

We define (instantaneous) projectors on to the spin-up 
and spin-down states as: 



f^: = i®\v+{t)){v+{t)\ 
n = i®\v_{t)){v.{t)\ ' 



(26) 



where / acts on the cantilever. The probability of spin- up 
is defined as 



and, the probability of spin-down is 



(27) 
(28) 
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The weight of the spin-up and spin-down wave packets 
is given by Vu and r^^ respectively. Though we can de- 
rive the value of rd from straightforwardly using equa- 
tion (|28|) . we use rd to represent spin-down probability 
for notational convenience. 

To keep track of the cantilever degree of freedom, we 
define variables that capture the mean position and mo- 
menta of the two wave packets as: 

:= {ZV^)/ru, Pu := {pVt)/^u ] 

. (29) 

Zd := {ZVi)/rd, Pd := {pVi)/rd ) 

We can now apply Eq. (|25]) to derive first-order mo- 
ment equations: 

dZu/dt = Pu/m, 

dpu/dt = -{muj^^2T)Zu + f(t), 

dZd/dt = Pd/m, 

dpd/dt = -{muj^-2r)Zd + f{t), 

dru/dt = 0. (30) 

The equations derived above close, which means that 
higher-order moment terms can be ignored. The ex- 
pected position of the cantilever is obtained by the iden- 
tity 

(Z) = TuZu ^ TdZd. (31) 

We observe that the moment equations actually are 
identical to the equations of motion of a driven (classi- 
cal) Harmonic oscillator; the resonant frequency of the 
two wave packets is shifted by an amount Ao; ^ F/mcj, 
from the natural resonant frequency uj of the cantilever. 
We also note that spin- up probability, Vu-, is a constant 
of motion. This is consistent with our derivation sec- 
tion mil Spin probabilities, however, are not a constant 
of motion for the non-unitary evolution considered in the 
next section. 



B. Open system evolution with continuous 
monitoring 

As discussed in Sec. HVl the cantilever-spin system is 
not isolated, and there are many sources of decoherence 
that affect both the cantilever and the spin. The sources 
include thermal noise, noise due to the process of contin- 
uous monitoring of the cantilever and noise in the mag- 
netic field. We model these decoherence effects in terms 
of Lindblad operators. 

There are three Lindblad operators that model the 
important sources of deocherence, two of which — Li 
and L2 — capture the decoherence acting on the can- 
tilever [Eq. mj)], while L3 models the noise on the spin 
[Eq. ([20]) ]. Li is the noise due to the thermal environ- 
ment and L2 is the noise due to continuous monitoring 
of the cantilever position. 



The evolution of the cantilever-spin system is given by 
a quantum stochastic master equation |13j]: 

dp = - (^-^ [Heff,p]dt 

^"^LkpLl - -{LlLk,p}dt 

+ V^d ((L2 - (L2)) p- p{Ll- (LI)) ) dWu 

(32) 

where Cd is the detector efficiency; dWt is a stochastic 
(white) noise process with the property: M[(iWt] = 
and dW^ = dt; and Heff is the effective Hamiltonian 
acting on the system [Eq. (p!8|) ]. 

As in the unitary case, we consider an arbitrary op- 
erator O acting on the system (with evolution described 
in Eq. ([32]) ). and derive its moment equations. Using 
Eq. ([32]) . the first-order moment equation for (O) = 
Tr{dp} is 

d{d) = '^{[Heff.d])dt 

+ ^(LtdL,)--({LtL,,0})^t 

fc=i 

+ ((OL2) - (6)(L2) 

+ (Lid) - (Lt)(6>) dWt. (33) 

We describe the (mean) position and momentum of the 
cantilever in terms of normalized wave packets with spin 
up or down, respectively [Eq. ([29]) ]. They are of the form: 

{0)„ = {dPt)/^« and (O), = (dn)/^d, (34) 

where Tu and Vd are the spin probabilities as defined in 
equations ([27]) and ([28]) . As in the unitary case, the mean 
(O) is defined as 

(O) = r„ (d)„ + u (0)rf. 

The system evolution described in Eq. ([32]) is stochas- 
tic. This means that the moments like {O)^ and {O)^ 
evolve stochastically as well; the moment equation ([33]) is 
a stochastic differential equation, unlike its unitary coun- 
terpart Eq. ([25]) . These equations follow stochastic cal- 
culus described by Ito rules [l5|. To give the evolution 
of (O)^ and (O)^, we apply Ito rules and get, 

d(0)„ = d(0)/r„ - {d)dr/rl - d{d)dr/rl 

+ {6){drf/rl (35) 

and similarly for (O)^. 
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1. First-order moment equations 

We now get to the task of deriving moment equations 
for the cantilever and the spin. As in the unitary case, 
we first derive the equation for spin-up probabihty = 
{V^)^ before we derive equations for the wave packets 
that describe the cantilever motion. 

Using Eq. ([33|) . we derive the equation for spin- up 
probability 



dvu = /^s (1 - 2ru)dt 



Zd)dWt. (36) 



Here, Zu and Zd are mean positions of the cantilever 
wave packets for spin up or down, respectively; tZs is 
the rate of spin noise [Eq. (f2Q|) ]: and A2 is part of the 
Lindblad operator for cantilever decoherence due to con- 
tinuous measurement [Eqs. ([16]) and (p!7|) ]. We used the 
identity in Eq. (|3T]) in our derivation. Note that unlike 
the unitary case, Eq. (|3Q|) . the spin- up probability Vu is 
not a constant of motion. However, since r^^ = (1 — r^^), 
the single variable is enough to capture the evolution 
of spin probabilities. 

The motion of the cantilever is described in terms of 
two wave packets corresponding to the spin being in the 
up or down state, respectively. Each wave packet is de- 
scribed by its position and momentum. To derive the 
moment equation for these variables, we apply Eqs. (j35j) 
and ([33]) . to obtain 



dZ„ 



and 



^{Pu/m) - 4edAl{AZ^)^rd{Zu 
- K.S {rd/vu) {Zu - Zd)dt 



Zd)] dt 



(37) 



dpu = {-{muj^ ^2T)Zu - 2-fPu + f{t))dt 
-Ks{rd/ru) {Pu - Pd)dt 
-4edAl{AR)^rd{Zu - Zd) dt 
^2^A2{AR)^dWt. (38) 

Similarly, 

dZd = [{Pd/m) + 4edAl{AZ^)^ru{Zu - Zd)) dt 
+ f^s iru/rd) {Zu - Zd) dt 
-^2^dA2{AZ^)ddWu (39) 

and 

dpd = {-{muj^ -2T)Zd - 2-fpd + f{t)) dt 
^ns{ru/rd) {Pu - Pd) dt 
^^edAl {AR)^ru{Zu - Zd) dt 
^2^A2{AR)^dWt. (40) 

In the unitary case, the first-order moment equations 
(|3Q|) formed a closed set. In the open system evolution 
however, the first order moment equations (|371 -HQ |) do 





/ 








~ ZuPu 




■■= 


-zj) 






- ZdPd^ 



not form a closed set, and depend on second-order terms 
like (AZ2)„ and {Ai?)„: 



(41) 
(42) 
(43) 
(44) 



where R was defined in (p!9]) above. Since the set of first- 
order equations do not close, we need to include second- 
order terms, and derive equations for them as well. 



2. Second-order moment equations and Gaussian 
approximation 

We will need equations for six second-order equations: 
(A^^)„, (Ap2)„, (Ai?)„, {AP)„ {Af), and (A^),. 
We begin with the derivation of (AZ^)^. 

To derive the moment equation, we first expand 
(AZ^)^ in terms of its definition (|4T]) . then apply Ito 
rules where necessary as shown in Eq. (|35]) for an arbi- 
trary operator. Finally, we expand using the moment 
equation ([33]) to get: 



d{AP)^ = {2{AR)Jm + Bin') dt 
-AeaAl ((AZ^)^ + (Z„ - Za) {AZ\) dt 
-^,(r,/r„) ((AZ2)„ - (AZ2), 

- {Zu - Zif ) dt 

+ 2^A2{AZ'')^dWt. (45) 

We have used some of the first order equations in the 
above derivation. 

We derive the equation for (AR)^, defined in Eq. 
by the same procedure used to derive Eq. ()35)) . and ob- 
tain: 

d{AR)^ = -{mJ +2r){AZ')^dt 
+ {{Ap\/m - 2 7{Ai?)„) dt 

-4e,Al ((Ai?)„(AZ2)„ 

+ rd(Z„-Zd){{AZ2,Ap}>„) dt 

-K,(rrf/r„) ((A^)„ - (Ai?), 

- {Zu - Zd){pu - Pd) ) dt 
+ y^A2 {{AZ^ Ap})„ dWt. (46) 



Note that a new second order term {Ap')^ appears in 
the above equation, and is defined as 



i^p\ ■■= {ip\ - pI) ■ 



(47) 
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Following the same procedure, we construct the moment 
equation for (Ap^)^: 



d{Af)^ = EUAlh'dt 
-4 ((mt^V2 + r) {AR)^+^{Ap\ 



dt 



■2edAi 



(2(Ai?): 



+ uiZ^-Zd) ({AZ,Ap2})^j dt 
■K«(rd/r„)((Ap2)„ - {Ap2)^ 

- {Pu - Vdf) dt 
-^dA2{{AZ,Af})^dWu (48) 



Similarly, we derive equations for the second-order mo- 
ment terms of the spin-down wave packet: 

d{AP)^ = {2{AR)Jm + B^h^'^ dt 

-4edAl ({AZ2)^ - r„(Z„ - Z^) {AZ^)^) dt 

+ irjra) ((Ai\ - (AZ^), 

+ {Zu - Zdf ) dt 



■2^A2{AZ^)^dWt, 



(49) 



d{AR)^ = - {moj^ -2T) {AP)^dt 
+ ({Ap')Jm - 2 7{Ai?>,) dt 

-4e,Ai {{AR),{AZ^), 

- ru{Zu - Zd){{AZ\Ap})^ dt 

+ «.(r„/rd) ((Ai?)„ - {Ai?>, 

{Zu - Zd) {pu - Pd)) dt 
+ ^dA2{{AZ\Ap})^dWu (50) 



d{Ap2)^ = Y^l^^Alh^dt 
-4 ((mu;V2 - T) {Ai?>, + 7(Ap'),) dt 
-2edAl (2{AR)l 

- r„ (Z„ - Zd) {{AZ, Ap2})^) dt 
+ Ks (r„/rd) ( (Ap2)^ - (Ap2)^ 

+ (pu - Pdf ) dt 

+ V^A2({Ai,Ap2})^dWt. (51) 

We have now derived equations for all the second order 
terms that appear in the equations. However, there are 
new terms, of third order, that appear in the second order 
equations (|^5l -l5T ]) : 

(AZ3)„_„ ({AZ, Ap2})„., and {{AZ\ Ap}>„_,. 



The second order equations (|45l -l5T ]) do not close as they 
depend on terms of third order. We now make the Gaus- 
sian approximation for the cantilever. Assuming that 
the wave packets were initially Gaussian, and that they 
remain approximately Gaussian at later times, then the 
third order terms vanish. The system of first-order equa- 
tions along with the second order equations now close, 
giving us a set of 11 coupled equations. 

We will now verify the validity of this approximation by 
numerically comparing the approximate solution to that 
of the full quantum stochastic master equation ([32]) . 



VII. NUMERICAL RESULTS 

We simulated the full quantum stochastic master equa- 
tion using a C + + quantum master equation library de- 
veloped by Brun and Shaw [l8|. We use a fifth order 
Runge-Kutta integrator [l^ to simulate the set of cou- 
pled equations (|36l -l4Q |) and (|45l -l5T ]) for the Gaussian 
approximation. 

We choose parameters based on those used by Berman 
et al. [9] and Brun and Goan [III, [HI- The values of 
parameters, in dimensionless units, are: 

h = uj = m = 1 
T] 0.6, e 100, 7 uj/Q 10-^ 0.85 
Al = 0.2, Bi = bx 10-^ A2 = 0.07, B2 = 0. 

(52) 

We assume the physical units in which uj = 10^ s~^, m = 
10-^2 kg, 7^ = 3 X 10^ T/m, and e = 300 /iT, consistent 
with current experiments. The values of Ai and Bi are 
different from the ones used in [11]. We assume that, in 
dimensionless units, 

ksT = 10^ 

in our present calculations. This value of ksT corre- 
sponds to a temperature of about 10 mK; though this is 
lower than the temperature accessible currently (which 
is around 300 mK [5]), given the considerable progress 
made on the experimental front in the last few years, 
the temperature we have used should be accessible in the 
near future. 

The output of the homodyne measurement that mea- 
sures the cantilever position is given by Eq. (|22]) , and can 
be expressed (in rescaled units) as 



Ic{t)dt = {Z)dt - A dWf 



where 



A 



50 



(53) 



(54) 



Here, Q is the quality factor of the cantilever, and the 
other parameters are due to the cantilever-light interac- 
tion, which is part of the optical interferometry process 



9 



(see Sec.|V]). Since the driving force apphed to the can- 
tilever is part of a feedback process, the driving frequency 
is resonant with the cantilever frequency. We know that, 
in the resonant case, the steady state amplitude of a Har- 
monic oscillator is proportional to its quality factor [l^ . 
In our simulation, we assume that the cantilever is in its 
steady state before it interacts with the spin, and that 
the steady state amplitude is, in dimensionless units, 50, 
which corresponds to 32 nm in physical units (in Q , the 
peak amplitude is 16 nm and the quality factor is 50, 000); 
due to this rescaling of amplitude, the cantilever quality 
factor Q appears in the above equation (|5^ . and we set 
P = -50/Q in Eq. ([22]). 

We can find the value of the coefficient A straightfor- 
wardly, by noting that (in dimensionless units) 



and 



Sf^E 



7c 



1.9 X 10^ 



1.4 X 10^ 



See [ll| (Sec. VII) for more information. 

The positive gain controlled feedback function f{t) is 
defined in Eq. (j23j), and depends on the feedback gain g 
and the pre-determined set point amplitude AMP. In 
our simulation, we set these parameters to 



g = 0.0001 and AMP = 50. 



(55) 



We chose a small value for g to avoid numerical instabil- 
ity when simulating the fully quantum stochastic master 
equation. However, the accuracy of the approximation is 
unaffected by larger values of g. 



exhibits this behavior because the value of His is relatively 
high (= 0.001), leading to multiple 'ffips' in spin proba- 
bility. 



B. OSCAR MRFM and single-spin measurement 

We now address the question of how high can be, 
for OSCAR MRFM to be a useful single-spin measure- 
ment device. Henceforth we use the numerically efficient 
Gaussian approximation equations to describe the OS- 
CAR MRFM system. 

In Figure [3l we plot the time evolution of spin- up prob- 
ability for different values of tZs- (0 10"^^ (^0 10~^ and 
10~^ in dimensionless units (16 Hz, 1.6 Hz and 160 mHz, 
respectively, in physical units); the parameters and the 
initial conditions are the same as in Fig. [2j We observe 
that the number of spin ffips decreases with the value of 
tvs' Intuitively, one can see from Fig. [3] that smaller val- 
ues of hv g are preferable, as we must track the cantilever 
position for a sufficient length of time to determine the 
orientation of the spin with required accuracy. If tZs is 
higher, the likelihood of a spin ffip increases, corrupting 
the samples used to determine the spin orientation. 

In the OSCAR protocol, we measure the spin by mea- 
suring the resonant frequency of the cantilever; more pre- 
cisely, the shift in the cantilever frequency. To determine 
the cantilever frequency numerically, we use the Fast 
Fourier Transform (FFT) algorithm dfourl from fv^ . 
The (numerical) frequency resolution of this FFT is 



A/: 



1 



TV At Tsampli 



(56) 



mg 



A. Accuracy of Gaussian approximation 

We now numerically demonstrate the accuracy of the 
Gaussian approximation by comparing it to a fully quan- 
tum stochastic master equation (|32]) . The parameters 
used to generate these plots are given in equations ([52]) 
and ([55]) : the spin noise rate, z^^, is set to 0.001. We as- 
sume that the cantilever and spin are initially decoupled; 
we also assume that the spin is in an equal superposition 
of up and down state, while the cantilever is at its lowest 
position {-AMP). 

In Figure [21 we plot the time evolution of the cantilever 
position and that of the spin-up probability; Figs. 2(a)[ 
2(b) and 2(c)[ shows the time evolution of the expected 
position of the cantilever (Z), as deffned in Eq. ([3T]) . at 
different times, while Fig. 2(d) plots the time evolution of 
the spin- up probability = {V^)- Figure [2] shows that 
the Gaussian approximation tracks the quantum stochas- 
tic master equation very closely. 

The spin probability does not relax completely (to ei- 
ther or 1); after initial ffuct nations, the spin probability 
gets closer to 1 (spin- up), the spin 'ffips' and the proba- 
bility is then closer to zero (spin-down), and ffnally, an- 
other 'ffip' takes it back closer to 1 (spin- up). The spin 



where N is the number of samples. At is the sample 
spacing, and A/ is the frequency resolution. TsampUng 
is the total sampling time. 

We assume that the sample spacing At is ffxed. The 
shift in the resonant frequency is proportional to F 
[Eq. ([T4]) ]. For the parameters chosen in our simulation 
[Eq. ([52|)], the value of F is 1.8 x 10"^. Thus, the fre- 
quency resolution to be achieved should satisfy 



A/ < r/27r ^ 2.8 X 10" 



(57) 



The sample spacing used in our simulation is At = 0.02; 
setting the frequency resolution A/ = 10~^, and using 
Eq. ([56]) . one can easily derive N = b x 10^. It is numer- 
ically convenient to the FFT algorithm if is a power 
of 2; we thus choose = 2^^ (the power of 2 closest to 
5 X 10^) for our simulation, and the frequency resolution 
becomes A/ = 9.5 x 10~^. 

Given the sample spacing. At, and the number of sam- 
ples, A^, the total sampling period in dimensionless (time) 
units is 



Ts. 



ampling 



NAt « 10,500. 



(58) 



The time scale of spin noise, k.^ ^, must be longer than the 
sampling period TsampUngi foi" OSCAR MRFM system 
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FIG. 2. (Color online.) Comparing full quantum SME and Gaussian approximation. We plot (Z) as a function of time for 
both the quantum SME as well and the Gaussian approximation. The parameter values used in this simulation are defined 
in equations (|52)) and (|55)) : the rate of spin noise Ks is set to 0.001. We assume that the cantilever-spin system is initially 
decoupled, with the cantilever at the lowest point in its oscillation, while the spin is in an equal superposition of up and down 
states. In Figs. 2(a) 2(b) and 2(c) we plot the expected position of the cantilever (Z) at different times, and in Fig. 2(d) we 
plot the evolution of spin-up probability = (^t)- 



to be useful in single-spin measurement. Thus, we set 
= 10~^ in our simulation. 

In Figure |4] we show the frequency shift achieved in 
the OSCAR protocol. We plot two different trajecto- 
ries: one in which the spin relaxes to the spin-up state, 
and another in which it relaxes to the spin-down state. 



Figure 4(a) plots the time evolution of the spin-u p prob- 
ability Tu for the two trajectories, and Fig. |4(b)] displays 
their corresponding Fourier amplitudes. It is clear from 
Fig. |4(b)| that there is a marked shift in resonant fre- 
quency: to the right in the case of spin-up, and to the 
left in the case of spin-down; we also observe that this 
shift is proportional to F (= 1.8 x 10~^). (Note that 
the natural resonant frequency of the cantilever, in di- 



mensionless units, is = 1.) A shift to the right in the 
resonant frequency of the cantilever implies that the spin 
is in the spin-up state, while a shift to the left indicates 
that the spin is in the spin-down state. 

The value of sampling period, TsampHng^ in Eq. (|58]) is 
in dimensionless units. In our calculation of parameter 
values for our simulations, we assumed that the cantilever 
frequency, in physical units, is 



phy 



phy 



/27r ^ 16 kHz. 



Using this frequency, the duration of sampling, in physi- 
cal units, is 

"^Sampling = ^Sampling / fc^^ ^ 656 mS, 



= 10" 
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FIG. 3. (Color online.) Spin probability for different Ks. We plot the time evolution of spin-up probability, ru — (V^), for 3 
different values of spin noise rate Ks- (i) 10~^, (ii) 10~^, and (m) 10~^. The number of spin flips decreases as Ks is reduced. 
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FIG. 4. (Color online.) Frequenc y shif t in OSCAR MRFM. We show the frequency shift in the OSCAR MRFM protocol, for 
two different trajectories. Figure 4(a) shows the time evolution of the spin-up probability, r^, for the two trajectories; in the 



first, the spin relaxes to its up state, while in the second, it relaxes to its down state. Figure 4(b) plots t he Fo urier amplitude 
(in arbitrary units) as a function of cantilever frequency, cj, corresponding to the two trajectories in Fig. 4(a) The spin noise 
rate k,s used in this plot is 10~^. The Fourier amplitude is calculated using a standard FFT algorithm; the number of samples 
is N = 2^^, and the sample spacing is At = 0.02. 
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and the shift in the cantilever frequency is 
YPhy = r X fP^y ^ 29 Hz. 

Thus, it takes about 650 ms for the OSCAR protocol to 
determine a shift of around 30 Hz in cantilever frequency, 
and consequently, to ascertain the orientation of the spin. 
The time scale of the spin noise must be longer than the 
sampling duration to use OSCAR MRFM as a single-spin 
measurement device. 



VIII. CONCLUSIONS 

In this work, we modeled the OSCAR MRFM opera- 
tion including decoherence effects on the cantilever (both 
thermal and monitoring), as well as spin noise due to 
magnetic sources. We describe the evolution of the sys- 
tem using a quantum stochastic master equation, the 
most general description in a Markovian framework. We 
then simplified this description using a series of approx- 
imations, and arrived at a semi-classical description of 
the system based on a Gaussian approximation of the 
cantilever state. 

We numerically compared the Gaussian approxima- 
tion to the fully quantum stochastic master equation, 
and found that the Gaussian approximation tracked the 
fully quantum stochastic master equation very closely. 
Thus, we conclude that the Gaussian approximation of 



the cantilever is useful for our range of parameter values, 
and that the MRFM system implementing the OSCAR 
protocol can be described by a closed set of 11 coupled 
stochastic differential equations. 

Finally, we used the Gaussian approximation equations 
to numerically verify OSCAR as a useful single-spin mea- 
surement protocol. The critical element in the OSCAR 
protocol is the time (sampling) it takes to determine the 
shift in resonant cantilever frequency. For the parame- 
ter values chosen in our simulations, we found that the 
sampHng duration is around 650 ms, and the frequency 
shift is about 30 Hz. The sampling duration thus places 
a bound on the time scale of spin noise that the OSCAR 
protocol, implementing a single-spin measurement, can 
tolerate. 
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